Statistics of small scale vortex filaments in turbulence 
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We study the statistical properties of coherent, small-scales, filamentary- like structures in Tur- 
bulence. In order to follow in time such complex spatial structures, we integrate Lagrangian and 
Eulerian measurements by seeding the flow with light particles. We show that light particles pref- 
erentially concentrate in small filamentary regions of high persistent vorticity (vortex filaments). 
We measure the fractal dimension of the attracting set and the probability that two particles do 
not separate for long time lapses. We fortify the signal-to-noise ratio by exploiting multi-particles 
correlations on the dynamics of bunches of particles. In doing that, we are able to give a first 
quantitative estimation of the vortex-filaments life-times, showing the presence of events as long as 
the integral correlation time. The same technique introduced here could be used in experiments as 
long as one is capable to track clouds of bubbles in turbulence for a relatively long period of time, at 
high Reynolds numbers; shading light on the dynamics of small-scale vorticity in realistic turbulent 
flows. 

PACS numbers: 47.27.-i, 47.10.-g, 47.ll.-j 



Vorticity dynamics, in general, and vortex filaments, 
in particular, have been the subject of many theoret- 
ical, phenomenological and experimental studies from 
the turbulent community, both applied and theoretical 
[Hi S S 13] ■ The transport of particulate by fluid flows 
is an ubiquitous phenomenon in nature and in indus- 
trial applications alike. According to their inertia prop- 
erties particles respond differently to fluctuations of the 
adverting -Eulerian- velocity field producing locally non 
homogeneous concentration, a phenomenon dubbed as 
preferential concentration 0, @| . As an example of such 
differential response, it has been clearly shown in the past 
that heavy particles tend to be expelled from high vor- 
ticity regions, while light particles tend to concentrate 
in regions where vorticity is higher Thanks to their 
strong property to concentrate in high vorticity regions, 
very light particles (i.e. small bubbles in water) have 
been used to visualize small scale vortex filaments 0, [3 . 
The strong tendency of light particles to concentrate in 
vortex filaments has been quantified also measuring the 
fractal dimensions of their dynamical attractor Q . Sim- 
ilar phenomena, based on complex response of micro- 
scopic hydrogen particles in quantum fluids, have also 
been exploited recently to visualize quantized vortices 

M- 

In the present work we will focus on the dynamics of 
light particles and will study in detail the connection be- 
tween their dynamics and the one of small scale vortex 
filaments 0. 

The data used for this study come from Direct Numer- 
ical Simulation (DNS) of 3D fully periodic Navier-Stokes 



eqs plus particles. Indeed, together with the Eulerian 
field we integrated the Lagrangian evolution of particles 
by mean of one of the simplest, yet nontrivial, model of 
dilute, passively advected, suspensions of spherical par- 
ticles as described in Refs. 



n|,ll2| : 
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(1) 



In the above equation x{t) and v{t) denote the parti- 
cle position and velocity respectively, r ~ a^/(3/3i^) is 
the particle response time, a is the particle radius and 
St = Tp/Tji is the Stokes number of the particle and 
/3 = 3pf/{pf + 2pp) is related to the contrast between 
the density of the particle, pp, and that of the fluid, pf. 
The incompressible fluid velocity u{x,t) evolves accord- 
ing to the Navier-Stokes equations : 



DtU = dfU + u ■ Vm = — Vp + v/S.u + / . 



(2) 



where p denotes the pressure and / an external forcing 
injecting energy at a rate e = (u ■ f). Eqn. ^ is evolved 
by means of a pseudo-spectral code for the fluid part with 
a second order Adams-Bashforth integrator, also used for 
the -millions of- particles evolving according to the dy- 
namics (H)), where the fluid velocity at the particle po- 
sition, u(x{t).^f)j wa s obtained by means of a tri-linear 



interpolation 



Energy was injected maintaining 



constant the spectral content of the first two shells in 
Fourier space. We report data coming from two sets of 
simulations with — 128^ and A^^ = 512^ collocation 
points, corresponding to Rca — 80 and Re^ — 180, re- 
spectively. We focus mainly on very light particles, in the 
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FIG. 1: PDF of halving-times (see text) for various density 
ratios, (3, at fixed St ^ 0.5, conditioned on the initial value 
of the vorticity magnitude to be larger than a given thresh- 
old. As P increases the PDFs show larger and larger tails, 
suggesting that lighter particles are probing regions of high 
and stable vorticity. 



limit of pp ^ (/3 — > 3 ) and on tracers evolving with 
the local Eulcrian velocity field x{t) = u{x{t),t). 

Inertial particles are not distributed homogeneously in 
the volume, centrifugal forces tends to concentrate light 
particle inside strong elliptical regions (with high vortic- 
ity) and heavy particles in hyperbolic regions, typical of 
intense shear. One thus expects different temporal corre- 
lations between particle trajectories and the underlying 
topology of the carrier flow. The local topology of the 
Eulerian flow is defined in terms of the symmetric and 
anti-symmetric component of the gradients Aij ~ diUj, 
hence there is an intimate link between the statistics of 
energy dissipation and/or enstrophy with particles evo- 
lution [l^, EEj [13 • To give an idea of this effect, we 
show in Fig. [T]the conditional probability distribution of 
the halving times, of the vorticity magnitude along 
particle trajectories. The vorticity halving time is de- 
fined, given a time t, as the first time lag after which 
the vorticity has becomes 1 /2 larger or smaller than the 
initial value, Ti/2(t) = min(T\uj{t + T)/uj{t) = 1 ± 1/2). 
The PDFs are conditioned in such a way that the vor- 
ticity magnitude at the reference time t is greater than 
a given threshold (for the case in Fig. [1] we chose this 
value to be bcurms', indeed the shape of the PDFs does 
not change significantly for higher values of the thresh- 
old). As one can see from the inset of Fig. [l] for a given 
Stokes, St 0.5, at changing the density contrast, f3, 
one moves towards higher and higher probability to ob- 
serve long halving times, i.e. light particles tend to live 
in regions of high and stable vorticity. 

Stable vortex structures, chaotically advected by tur- 
bulent flows, should therefore play the role of attrac- 
tive sinks in the dynamics. Larger their lifetime will be. 



larger inhomogeneous bubble distribution will develop. 
To quantify better the clustering properties of light par- 
ticles, we start from studying the statistics of the largest 
Lyapunov exponent, Ai. It is known that clustering and 
inertia may affect the whole distribution of Lyapunov 
exponents. One would expects that inertia reduces the 
tendency for particles pairs to separate. Indeed, a small 
tendency toward a reduction of the Lyapunov exponent 
at increasing the Stokes number has been reported for 
heavy particles 18|, similar results have also been ob- 
tained in random flows [lij . Even more interesting is the 
study of the whole probability distribution of the largest 
finite-time Lyapunov exponent (FTLE), 71 (T) defined as: 
71 (T) = ^\og{R{T,t)/R{0,t)) where R{T,t) is the sep- 
aration of two particles at time t + T starting from a 
separation i?(0, t) at time t. For large times, the distri- 
bution of FTLE is expected to follow a large-deviation 
form: prili) exp[— rS'(7i)], where the Cramer func- 
tion, >5'(7i) is a non- negative convex function vanishing 
for 71 = Ai. In Fig. [2] we report the Cramer function at 
T = 19Qtjj for the case of fight particles f3 ^ 3, St ^ 1.2 
compared with the one obtained for tracers. As one can 
see there are two remarkable effects. First, the minimum 
for the case of light particles is achieved for a value much 
smaller than the one for the tracers, precisely ^it^j ^ 0.04 
for bubbles and 71 Tj, ~ 0.14 for tracers. Second light par- 
ticles show a remarkably high probability to have pairs 
that do not separate at all, even for long time lapses, i.e. 
there are many events in the left tail of the Cramer func- 
tion which have a negative FTLE. The global average 
properties, however, remains chaotic, i.e. with probabil- 
ity one all couples separate at long times. The Kaplan- 
Yorke dimension for the family of light particles shown in 
Fig. [2]is Dky = 1.3±0.3 0. The observation that such 
small-scales strong clustering properties are correlated to 
the topology of the fiow structures at the same scales 
suggests the possibility to use light particles to study 
statistical properties of small-scales vortex filaments in 
turbulence. 

The identification of small scale vortex filaments is an 
extremely daunting task. An analysis of the Eulerian 
fields would require the detection of isosurfaces of vor- 
ticity (larger than some prescribed threshold): a method 
problematic if not for the larger and most intense struc- 
tures . Furthermore to analyze the temporal evolution 
of vortex filaments one would also need to track three di- 
mensional structures not only in space, but also to follow 
their evolution in time by repeating the same analysis at 
each Eulerian time step. Our original proposal here is to 
use multi-particles correlation to extraordinary enhance 
the signal-to-noise ratio associated to the identification 
of vortex filaments. In Fig. [3]we show the trajectories of 
several light particles which are attracted, in presence of 
favorable pressure gradients, into a vortex filament and 
that then separate again once the vortex filament disap- 
pears. Only the favorable pressure gradient, which tends 
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FIG. 2: Cramer entropy for light particles and tracers. The 
function, for the case of light particles, is shifted towards lower 
values (the maximum Lyapunov exponent where 5(71) attains 
its minimum is smaller than the one for tracers) and shows 
larger tails (thus suggesting a higher degree of intermittency) . 
There is a non-negligible fraction of events (left tail of the 
red curve) with couples of particles that do not separate at 
all (absent in the tracers case). 
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FIG. 3: Evolution of the trajectory of different light particles 
entering in a small scale vortex filament and then separating 
again. In the inset it is shown the behavior of the momentum 
of inertia of the particles cloud. At entering (exiting) from 
the vortex filament an abrupt -notice the log scale- decrease 
(increase) of the signal allows to sharply define the life-time 
of the vortex filament. 



to concentrate bubbles in the vortex core, permits such a 
strong clustering. Once the vortex filament breaks down, 
particles are released and separate explosively. 

One can define the signal-to- noise (s/n) ratio for an 
event with n particles in a volume where on average there 
are n as: SNR = n/n. Our goal here is to identify an 
observable which is very sensitive to the presence of small 
scale vortex filaments. With such an observable it will be 
possible to fix a threshold to define a birth and death time 
for the vortex filaments. We now outline our procedure 



in detail. We take a snapshot of light particles configu- 
ration (i.e. positions) at a time roughly in the middle of 
our numerical simulation. We divide our simulation vol- 
ume into small cubes, say of size 4'^ Kolmogorov scales, rj, 
and we look for those volumes with larger particle counts. 
The particles residing in a volume will form what we call 
a bunch. We then consider the full trajectory of the par- 
ticles within the several bunches that we identified. Obvi- 
ously the number of particles in the different bunches will 
not be the same. We identify M bunches with the i-th 
bunch formed by Ni particles, then we define the center 
of mass of the bunch i as a;J.,„(t) ~ Sj^i '^j(^) 

and the momentum of inertia for the same bunch of par- 
ticles as: Mi = Tfjli i^jit) - Kmi^^)) ■ The physical 
interpretation is clear, the smaller the momentum of iner- 
tia of the bunch, the closer the particles. The important 
observation here is that this quantity is very sensitive 
to vortex filaments and displays an extraordinarily high 
signal/noise ratio (see inset of Figure [S]). To understand 
this point one has to consider the fact that we find eas- 
ily hundreds of particles at distances smaller than O.I77, 
while for an uniform distribution one would expect to 
find 1.4 • 10"'^ particles per 77'^. The probability to find a 
finite number of particles (even if just a few, say 3-^4) 
close inside 77'^ is so small that if this happens it is al- 
most surely associated to the presence of confining forces 
keeping the particles close-by (the probability could be 
estimate by means of the Poissonian distribution). 
Another remarkable feature that makes the momentum 
of inertia of the bunch an extremely useful quantity to 
identify vortex filaments is the rapidity with which light 
particles in the neighborhood of a forming vortex fila- 
ment arc attracted into it. In the inset of Fig. [3] one can 
indeed see that particles initially separated move closer, 
remain very close to each other for some time, and then 
separate again. The extreme sharpness (rapidity) with 
which the momentum of inertia drops and then raises 
again allows to define the life-time of the vortex filament 
in a robust way. The life-time estimates are not much 
sensitive to the chosen threshold (in Fig. [3] the threshold 
has been set to the value 1, changing the value of the 
threshold is accounted into the error bars e.g. in Fig. [4]). 
By analyzing the statistics of all bunches, we can make 
an histogram of vortex filaments life-times, allowing for 
the first time to assess, in a quantitative way, the sta- 
tistical properties of these extreme events. In Fig. 2] we 
show the pdf of vortex filaments life-times for Rex = 180 
and 80 which happens to be an exponential with a de- 
cay rate which can be estimated of the order of 25r^ and 
17r^, respectively. As already noticed in Q we do ob- 
serve events as long as the integral time Tl, that in our 
simulation was estimated to be of the order of 5Gr,j. An- 
other interesting question is how particles bunches are 
formed, if this happens by means of a sudden attraction 
of close-by particles, or as the result of a sequence of suc- 
cessive -and rapid- bunch merge. The same question can 
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FIG. 4: Probability distribution function of the life-time of 
vortex filaments. PDFs at two different Reynolds number, 
R\ — 80 and R\ — 180 respectively: in both cases it turns 
out that the PDF is fairly well fitted by an exponential (with 
a decay factor, for the highest Rx, of around 20t,,). In the 
far tail of the PDFs there are events whose life-time is of 
the order of the integral time. Error bars reflect the effect of 
changing the threshold used to determine the vortex filament 
time of life (see text for more details) . 



be posed regarding the "decay" of a vortex filament. How 
do particles separate? To answer this question one can 
introduce the following quantity measuring the average 
distance between two particles in a bunch: 



D 



2 ^' 



(x.it) - xk{t)y 



(3) 



Should all particles be close then the histogram of the 
distances would be centered around a small unique value 
with a variance connected with the bunch size. In 
the event of a bunch being splitted into two separating 
(smaller) bunches, one would expect to see the appear- 
ance of a peak at some finite distance and a shift of the 
position of this peak with time. This is exactly the case 
in Fig. [5] where the pdf of for a given bunch is shown 
for 4 consecutive instants in time. A clearly visible peak 
can be identified and its position is shifting towards larger 
distances. 

One of the most intriguing features of fluid dynamics 
turbulence is the presence of long living coherent struc- 
tures at small scales. The quantification of the statistical 
properties of these Eulerian structures has always proved 
to be one of the most difficult statistical analysis in tur- 
bulence due to the extremely low signal to noise ratio and 
to the need of following the evolution of structures for as 
long as a large scale eddy turnover time. We showed, for 
the first time, that the coherent dynamical properties of 
clouds of very light particles (bubbles) can be used to in- 
troduce an observable extremely sensitive to small scale 
vorticity filaments. We used this observable to quantita- 



FIG. 5: The probability distribution function of the aver- 
age distance (squared) between pairs of particles in a bunch 
clearly showing a fragmentation process: a vortex filament is 
splitted in two and the bunch of particles is splitted in two 
bunches separating in time (the distance between the bunches 
is associated with the moving position of the peak) . 



tively measure the probability distribution of vortex fila- 
ments life-times. We also show that the decay of a vortex 
filament into smaller chunks is the way particles get re- 
leased from within such structures. The same observable 
introduced here could be used in experiments as long as 
one is capable to track clouds of bubbles in turbulence for 
a relatively long period of time, at high Reynolds num- 
bers. Vortex dynamics is of particular importance also to 
study possible blow-up of Euler and Navier-Stokes equa- 
tions 



22l | as well as super-fiuid dynamics 2l[. Similar 
ideas could abridged to other domain to fortify statisti- 
cal signals. 
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